VT asymmetry in viscous fluids with balanced inflow and outflow 
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In recent years, open systems with equal loss and gain have been investigated via their symmetry 
properties under combined parity and time-reversal iVT) operations. We numerically investigate 
T^T-symmetry properties of an incompressible, viscous fluid with "balanced" inflow-outflow config- 
urations. We define configuration-dependent asymmetries in velocity, kinetic energy density, and 
vorticity fields, and find that all asymmetries scale quadrat ically with the Reynolds number. Our 
proposed configurations have asymmetries that are orders of magnitude smaller than the asym- 
metries that occur in traditional configurations at low Reynolds numbers. Our results show that 
T^T-symmetric fluid flow configurations, which are defined here for the first time, offer a hitherto 
unexplored avenue to tune fluid flow properties. 



Introduction. In a seminal paper [1 Carl Bender and 
co-authors showed that all eigenvalues of a class of non- 
Hermitian Hamiltonians are purely real and bounded 
below. For a Hamiltonian Hr^ = —d^ + x^(zx)^, with 
complex potential Vrj{x) = x^(ix)^, they showed that 
when ?7 > 0, the eigenvalue spectrum is purely real 
whereas for 77 < 0, it has a finite number of purely 
real eigenvalues and the rest occur in complex conju- 
gate pairs. Since then, the properties of Hamiltonians 
that are invariant under the combined operation of par- 
ity and time-reversal {VT) have been extensively stud- 
ied in the continuum [2]-[4] and lattice [5-10 models. 
For a T^T- symmetric Hamiltonian, the region of param- 
eter space where all eigenvalues e^ are real, and the 
eigenfunctions '^a(^) are simultaneous eigenfunctions of 
the combined T^T-operation, is called the T^T-symmetric 
phase; the emergence of complex conjugate eigenvalues 
that occur when the system leaves the T^T-symmetric 
phase is known as T^T-symmetry breaking. This tran- 
sition is characterized by the asymmetry parameter de- 
fined as ppt = ^x J dx\iljx{x) — '^^(— x)|. In the VT- 
symmetric phase pPT = 0, whereas in the 'PT-broken 
phase, this asymmetry parameter is nonzero. The con- 
verse of this observation, proven recently, shows that if 
a non-Hermitian Hamiltonian has purely real spectrum, 
it must be T^T-symmetric with an appropriately defined 
parity operator that satisfies V^ = 1 pTl . 

Recent experiments have shown that, far from being 
a mathematical curiosity, T^T-symmetric Hamiltonians 
naturally represent open, non-equilibrium systems. VT 
symmetry breaking has been observed in evanescently 
coupled optical- waveguides, where light is absorbed in 
one and amplified in another p!2Hl5] . coupled electrical 
circuits [TB^, and coupled pendulums ^. These studies 
have shown that open systems with spatially separated, 
balanced sources and sinks undergo a transition from zero 
asymmetry state to a state with ppT 7^ [181 E] • 



The notions "source-sink", "gain- loss", and "inflow- 
outflow" occur most naturally in fluid systems. However, 
symmetry properties of viscous flow are not straightfor- 
ward because the non-linear nature of underlying dynam- 
ics implies that superposition principle cannot be used 
to construct solutions with specific symmetries. Tradi- 
tionally, viscous flows are driven by upwind flow, pres- 
sure difference, or boundary movement and therefore the 
(steady-state) velocity profiles at the inlet and the outlet 
are, in general, not identical. In particular, flow symme- 
try properties in a system with identical inflow and out- 
flow velocity profiles remain unexplored [20 . We note 
that in the experimentally investigated T^T-symmetric 
systems [T3'-'T6 , the gain and loss were confined to elec- 
tromagnetic fields or energy [17 ; for a viscous fluid sys- 
tem, one must distinguish between the mass flow, defined 
by T^T- symmetric boundary conditions, and the energy 
flow properties, characterized by the fluid viscosity. 

In this paper, we introduce incompressible, viscous 
fluid configurations with T^T-symmetric boundary condi- 
tions and numerically investigate the symmetry proper- 
ties of resultant steady-state flow using lattice Boltzmann 
method. We define and calculate the T^T- asymmetries in 
velocity, kinetic energy density, and vorticity as a func- 
tion of Reynolds number and distinct configurations of 
"balanced" inflow and outflow. Our salient results are as 
follows: i) The asymmetries in balanced configurations 
are suppressed by orders of magnitude when compared 
with those in traditional configurations, ii) These asym- 
metries scale quadrat ically with the Reynolds number, 
iii) These results are valid across different balanced in- 
flow and outflow configurations. We emphasize that in 
all configurations - VT symmetric or traditional - the 
total mass flux at the inlet and the outlet are identical. 
Thus, the phrase "balanced inflow and outflow" implies, 
in addition, symmetry properties of velocity profile at the 
inlet and the outlet, as we define below. 
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FIG. 1. Two distinct balanced inflows and outflows possible 
in a square geometry. Panel (a) shows a T^T-invariant con- 
figuration where parity Vo corresponds to reflection across 
the square diagonal. Panel (b) shows a T^T-invariant config- 
uration where parity Ve corresponds to reflection through 
the origin. We use a parabolic velocity profile Ux{y) = 
Up[l — {2y/wY]0(w — 2\y\) at the inlet and outlet, where 0{y) 
is the Heaviside function. 

Model and Formalism. We start with the definition 
of T^T-symmetry for fluids, motivated by the notion of 
balanced inflow and outflow. Consider two square ge- 
ometries of width W, with inflows /out flows of width w 
and origin of the co-ordinate system at the center of 
the square. Both panels, (a) and (b), in Fig. IT] show 
"balanced" situations, but correspond to different defini- 
tions of parity operator. The left-hand panel (a) shows 
boundary velocity profile that is invariant under VT- 
transformations where the parity operator is defined as 
Vo{x^y) = {—y^—x). This (odd) parity operation, a 
reflection in the southeast-to-northwest diagonal, satis- 
fies detT^o = — 1 and Vq = 1. Under this transfor- 
mation, the velocity field inside the square changes as 

The right-hand panel (b) shows boundary velocity 
profiles that are invariant under T^T-transformations 
where the even "parity" operator is given by VE{x^y) = 
(—X, —y). In contrast to the standard definition of par- 
ity in even dimensions, detT^^; = 1 since it is just a 
two-dimensional 7r-rotation. The velocity field remains 
unchanged under this transformation, VET(ux^Uy) = 
(ux^Uy). Note that one-dimensional inflow-outflows are 
special cases of this geometry with a = (vertical) and 
a~^ = (horizontal). We point out that apart from the 
inflow and outflow locations, fluid velocities at all other 
points on the boundary vanish due to no-slip condition 
required by a viscous fluid, and that for a = {0, 1, oo}, 
the boundary conditions are T^T-symmetric with respect 
to both parity operators Vo and Ve- 

We define the dimensionless VT asymmetry in the 
steady-state velocity field u(r) for panel (a) configura- 
tion as 
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]^ j dv [m{Vov) - VoTu{rf , (1) 



and corresponding definitions for asymmetry pQ in the 
kinetic energy density T(r) = pu(r)^/2 where p is the 



fluid density, and the asymmetry pQ for the pseudo-scalar 
vorticity field 0Jz{'^) = V x u(r). Eq.Q is applicable for 
a < 1 whereas when a > 1, the prefactor in front of the 
integral changes to [2W'^ {aup)'^]~^ so that the equiva- 
lence between a ^ a~^ and the exchange of axes, x ^ y^ 
is preserved. For configurations in panel (b), correspond- 
ing asymmetries p^ are defined using the even "parity" 
operator Ve instead of the odd parity operator Vo • Note 
that, by construction, the boundary contribution to the 
asymmetry in all "balanced" inflow and outflow scenarios 
is zero. 

The fluid flow is governed by the Navier- Stokes equa- 
tion 



dtVi + (u • V)u = -Vp/p + i^V^u, 



(2) 



where v is the kinematic viscosity, and u(r, t) and p(r, t) 
are the velocity and pressure fields respectively. The 
Reynolds number is defined as Re = lipi^/i^ for a < 1 
(Re = aupw/v for a > 1). We use kinetic-based lattice 
Boltzmann method (LBM) to solve Eq.([2| through a val- 
idated code. Originated from lattice gas automata, the 
LBM has developed into a promising numerical scheme 
for modeling viscous fluids, with physical and computa- 
tion advantages to deal with flow complexity ^TJ [22] , 
The fundamental idea of the LBM is to construct sim- 
plified kinetic models for spatially and temporally dis- 
cretized distribution functions. The evolution of these 
functions incorporates the physics of mesoscopic pro- 
cesses so that the coarse-grained observables obey de- 
sired macroscopic equations [23 . The lattice Boltzmann 
equation reads 

fp{r + ep5ut + 5t) = Mr,t)-[Mr,t) - fp{r)\ /r (3) 

where //^(r,^) (/3 = 0, • • • , 8) are the single-particle dis- 
tribution functions with nine discrete velocities e/3 = 
{0, ±x, ±^, ±x ± y} used in the two-dimensional 9-bit 
square lattice model [23], /«^(r) are the equilibrium dis- 
tribution functions, r is the relaxation time determined 
by molecular collisions, and in turn, it determines the 
viscosity of the fluid. The hydrodynamic velocity field is 
given by 



u(r) = ^E^/^/z^w = E^/^//^w/E//^w- 



(4) 
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We use bounce-back boundary condition at the walls to 
ensure no-slip, and generalized bounce-back boundary 
condition [24] at the inlet and the outlet with a parabolic 
velocity profile having a centerline velocity Up or aup. 
Numerical Results. Let us start with the steady-state ve- 
locity profiles obtained for different configurations to get 
a feel for typical velocity fields (left-hand side in Fig. [2]) . 
Panel (a) has a = and shows a small velocity asym- 
metry as can be visually checked by comparing velocities 
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FIG. 2. Left figure: Typical steady-state velocity fields u(r) for different "balanced" geometries with w/W = 0.1; a small 
Re = 0.01; and a = (a), a = 1 (b), a = 0.58 (c) and (d). Note that due to the geometry, the velocity profile in (a) satisfies 
Ux{—x,y) = —Ux{x,y) and Uy{—x,y) = Uy{x,y), and in (c) satisfies Ux{y,x) = Uy{x,y) and Uy{y,x) = Ux{x,y) exactly. Right 
figure: For configuration in panel (a), dependence of the steady-state u(r) as a function of Reynold number. For Re = 0.01 
(A), the velocity profile is essentially symmetrical; for Re = 10 (B) and Re = 30 (C), the asymmetry is visible. When Re = 100 
(D) the asymmetry is accentuated by the emergence of vortices near the infiow region that are absent near the outfiow region. 



at points ±r. When a = 1, panel (b), the resultant ve- 
locity profile is reflection-symmetric about southwest-to- 
northeast diagonal, shown by the white dashed line, but 
has a weak asymmetry about the other diagonal, shown 
by white dash-dotted line. Panel (c) with a = 0.58 has 
a weak asymmetry across the diagonal denoted by white 
dash-dotted line and Panel (d) shows the complementary 
conflguration. Note that for panels (a), (b), and (d), the 
asymmetry pE associated with reflection through the ori- 
gin, shown by the central white dot, is relevant; for panels 
(b) and (c), the asymmetry po associated with reflection 
along the diagonal, shown by white dash-dotted line, is 
relevant. 

The right-hand side in Fig. [2] shows the steady-state 
velocity fleld for a vertical flow {a = 0) as a function 
of the centerline speed Up, or equivalently, the Reynolds 
number Re. As Re changes over four decades, the asym- 
metry in the velocity fleld at points ±r also increases. 
When Re = 0.01 (A), the even-parity asymmetry is 
barely noticeable; for Re = 10 (B) and Re = 30 (C), 
it is clearly visible; for Re = 100 (D), the asymmetry 
is accentuated by the emergence of vortices near the in- 
flow region that are absent near the outflow region. Note 
that the velocity fleld satisfles Ux{—x^ y) = —Ux{x^ y) and 
Uy{—x^y) = Uy{x,y) exactly, and thus guarantees that 
the net vorticity of the velocity fleld vanishes. 

Next, we consider the quantitative scaling of the rel- 
evant asymmetries with system parameters. The left- 
hand side of Fig. [3] shows the typical scaling of even- 
parity asymmetries p£;(Re,Qf) with Reynolds number Re 



for different values of a < 1. Panel (A) shows that the 
velocity asymmetry scales quadratically with the Reynolds 
number over four decades^ p'% = A'^{a)Re^; panels (B) 
and (C) show that the asymmetries in kinetic energy 
density and vorticity also scale quadratically with the 
Reynolds number, p^ = A^{a)Re^ and p^ = A^{a)Re^. 
Panel (D) shows that the prefactor functions for the 
even-asymmetry A* (a) increase monotonically with a for 
< a < 1. These results are valid for different inlet sizes 
or boundary velocity proflles. 

The right-hand side of Fig. |3] compares the asym- 
metries in balanced-outflow model (BO) with those of 
a traditional model with fully-developed-outflow (FDO) 
at the outlet. Panels (a) and (b) show the speed pro- 
flles at Re = 0.01 and w/W = 0.1 for the two cases 
respectively. Panel (c) shows the clear difference be- 
tween the steady-state velocity proflle Ux{x^y) at the 
outlet for the balanced-outflow case (solid line) and the 
fully-developed-outflow case (dashed line). These differ- 
ences between the PT-symmetric and traditional models 
are quantifled in panel (d). It shows that in the VT- 
symmetric case, as before, the relevant odd asymmetries 
scale quadratically with the Reynolds number (lines with 
open symbols). In a sharp contrast, however, asymme- 
tries in traditional models with fully- devloped- outflow are 
orders of magnitude larger and, over three decades, are 
mostly insensitive to the Reynolds number (lines with- 
out symbols). These results show that T^T-symmetric 
inflow-outflow conflgurations strongly suppress observ- 
able asymmetries compared to their traditional counter- 
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FIG. 3. Left figure: For panel (d) of Fig. [2] tiie dependence of even-parity asymmetries in velocity (A), kinetic energy density (B), 
and the vorticity (C) as a function of a shows that they scale quadratically with the Reynolds number, p^(Re, a) = A^{a)Ke^ 
for i = u,k,uj. Panel (D) shows the behavior of A^(a) for < a < 1. Right figure: Comparisons between balanced-outflow 
(panel a) and fully-developed-outfiow (panel b) speed profiles at Re = 0.01 and inlet size w/W = 0.1. Panel (c) shows the 
clear difference between velocity profiles Ux{x = VK/2, y) at the outlet for BO boundary condition (solid line; panel a) and FDO 
boundary condition (dashed line; panel b). Panel (d) shows that po(Re) for the BO case (lines with symbols) scale as Re^, 
whereas those for the FDO case (lines without symbols) are orders of magnitude larger. 



parts. 

Discussion. In this paper, we have developed the formal- 
ism of T^T-symmetric inflow-outflow configurations for 
incompressible viscous fluids. Motivated by the loss of 
reciprocity, pPT 7^ 0, as the signature of T^T-symmetry 
breaking, we have defined configuration-dependent asym- 
metry functions p(Re, a) for fluid velocity profile, ki- 
netic energy density, and vorticity. Although the non- 
linearities due to convective acceleration, the (u • V)u 
term in the Navier-Stokes equation, make it difficult 
to predict the symmetry properties of these variables, 
we have found that all asymmetries scale quadratically 
with the Reynolds number over a wide range of geome- 
tries and boundary velocity profiles. We emphasize here 
that p^(Re) and p^(Re) involve the square of velocity, 
whereas the kinetic energy density involves its fourth 
power. Thus, their identical quadratic scaling is can- 
not be ascribed to the functional form of the asymmetry. 
Similarly, the virtual independent of these asymmetries 
on the Reynolds number in the fully-developed-outflow 
case implies that this effect cannot be purely ascribed to 
the convective acceleration term. 

Our results show that solutions of incompressible, 
viscous flow problems with T^T-symmetric boundary 
conditions are highly symmetrical, particularly at low 
Reynolds numbers Re < 1, compared to their counter- 
parts with traditional boundary conditions. This in- 
creased symmetry of the steady-state solutions implies 
that such configurations are, generically, robust against 
the formation of vortices than the traditional configura- 



tions. Our results thus lay the foundation for the study 
of T^T-symmetric fluid systems. Note that the lack of a 
finite, nonzero threshold for the asymmetry functions is 
due to the purely dissipative, viscous term and nonlinear- 
ities in Eq.([2]); this is in contrast with the experiments 
on T^T-symmetric systems [TSHIZ] where both dissipa- 
tion and amplification are present. Formally, an exactly 
T^T-symmetric fluid system requires an antisymmetric 
viscosity v{Vy) = — i^(r). Since this is physically im- 
possible, one can consider a system where viscosity has 
an antisymmetric component u{x^y) = i^o + sign(x)Ai^. 
This location-dependent enhancement or suppression of 
viscosity is possible in ferrofluids [25l|26]. Thus, investi- 
gation of such a system may provide an example of truly 
'PT-symmetric fluid system. We also note that viscous 
fluid flow with porous walls has been extensively studied, 
although such problems do not have "balanced inflow and 
outflow" boundary conditions [27H29] . 

This work raises a number of questions: Does the 
quadratic scaling of asymmetry persist when the flow be- 
come transitional or turbulent at higher Reynolds num- 
bers? How does the asymmetry scale in a rectangular ge- 
ometry? Does the onset of turbulence occur at the same 
Reynolds number for a "balanced" geometry as it does 
for fully-developed-flow boundary conditions? Answers 
to these questions may inspire innovative flow control 
techniques, with implications to a wide variety of fields. 
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